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We study gravity coupled to a scalar field in spherical symmetry using loop quantum gravity 
techniques. Since this model has local degrees of freedom, one has to face "the problem of dynamics" , 
that is, difi^eomorphism and Hamiltonian constraints that do not form a Lie algebra. We tackle 
the problem using the "uniform discretization" technique. We study the expectation value of the 
master constraint and argue that among the states that minimize the master constraint is one that 
incorporates the usual Fock vacuum for the matter content of the theory. 



I. INTRODUCTION 



o 
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\ Loop quantum gravity is being explored in model situations of increasing complexity. There has been steady 

(~| ■ advance in treating homogeneous cosmologies [l|, an area of activity that has come to be known as loop quantum 

^ cosmology. There has also been progress in spherical symmetry in vacuum [3]. However, in all these cases one did 
not have to face the "problem of dynamics" , i.e. dealing with the non-Lie algebra of constraints of general relativity. 

' In homogeneous cosmologies there is only one constraint and it therefore has a trivial algebra. In spherical symmetry 

, special gauges were chosen that resulted in an Abelian algebra. In this paper we would like to study spherically 

^ ■ symmetric gravity coupled to a spherically symmetric scalar field using loop quantum gravity techniques. It is not 
known in this situation how to formulate the problem in a way that one ends up with a Lie algebra of constraints. A 

' total gauge fixing was introduced by Unruh but it leads to a non-local expression for the Hamiltonian. Here we 

(2jr)| will fix partially the gauge to eliminate the diffeomorphism constraint in order to simplify things. This still leads to 
'— a Hamiltonian constraint that has a non-Lie Poisson bracket with itself, involving structure functions. To treat this 

^ ' problem we will use the "uniform discretization" technique ^] . We will introduce a variational technique adapted to 

^ ■ the minimization of the master constraint (in the context of uniform discretizations one should probably refer to it as 

■ "master operator" since it only vanishes in the continuum limit). In the case that zero is in the kernel of the master 
constraint the technique yields the correct physical state in model situations. 

C The inclusion of scalar fields in spherical symmetry opens a rich set of possibilities to be studied including the 

'""^ formation of black holes, critical collapse, the emergence of Hawking radiation, among others. Here we will have much 

\^ more modest goals: to see how the complete theory approximates the vacuum state of the scalar field living on a fiat 

' space-time. An outstanding problem in a full quantum gravity treatment involving matter fields is the emergence of 

0^ , a vacuum state for the fields and what relation it may have to the ordinary Fock vacuum of quantum field theory in 

' curved space-time. We will apply the variational technique in the case of spherically symmetric gravity coupled to a 

J> , scalar field and show that it yields a vacuum state that is closely related to the Fock one. 

The organization of this paper is as follows: in section II we review the classical theory. In section HI we discuss 

. , the quantization of a spherical scalar field in a classical fiat space-time in order to have something to compare with 

■ the full case. In section IV we study the full quantization of gravity and the scalar field, using a variational technique 
" " ' to minimize the master constraint. We end with a discussion. 
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II. SPHERICALLY SYMMETRIC GRAVITY WITH A SCALAR FIELD: THE CLASSICAL THEORY 

Spherically symmetric gravity with the Ashtekar new variables has been studied in detail in Q and Here we 
present only a brief summary. One assumes that the topology of the spatial manifold is of the form E = i?"*" x S"^ . We 
will choose a radial coordinate x and study the theory in the range [0, oo]. The invariant connection can be written 
as, 

A = A^{x)A3dx + {Ai{x)Ai + A2{x)A2) de (1) 
+ {{Ai {x)A2 - A2 ix)Ai) sin 9 + A3 cos 9) dip, 

where Ax,Ai and A2 are real arbitrary functions on the A/ are generators of sm(2), for instance A/ = —iai/2 
where cr/ are the Pauli matrices or rigid rotations thereof. The invariant triad takes the form, 

E = E^x)A3sm9— + {E\x)Ai+E\x)A2)sm9— 
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+ {E\x)K2-E\x)K,)-^, (2) 

where again, , and E"^ are functions on 

As discussed in our recent paperQ and originally by Bojowald and SwiderskiQ, it is best to make several changes 
of variables to simplify things and improve asymptotic behaviors. It is also useful to gauge fix the diffcomorphism 
constraint to simplify the model as much as possible. It would be too lengthy and not particularly useful to go 
through all the steps here. It suffices to notice that one is left with two pairs of canonical variables E'^ , and E'^, 
and that they are related to the traditional canonical variables in spherical symmetry ds^ = A?dx^ + R^dil^ by 
A = E'^/y/\E^\, Pa = -^/\E^\K^, R = and Pr = -2y^\E^\Kx - E'^K^/^/\E^\ where Fa is the momentum 

canonically conjugate to A. 

In terms of these variables the diffcomorphism and Hamiltonian constraints for gravity minimally coupled to a 
massless scalar field are 0, 

Cr = {\E-\yKx-E'^{Kj -P^^' (3) 
1 



^ G 



2y/\W\ 2y/\W\ 8y/\W\EV 



2{EvY 2EV 



2^/\EF\Ef 'i'E'c 



(4) 



and since the variables are gauge invariant there is no Gauss law. We have taken the Immirzi parameter equal to one. 
We now proceed to partially fix the gauge by choosing E^ = x'^ {R = x in terms of the metric variables). One can 
solve the diffcomorphism constraint for Kx, 

= ,5) 

which yields the Hamiltonian constraint for the partially gauge fixed model as, 
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E-^ E'fKl 3a; x'^iE'^J ^ 
2^ " + 2^ ~ {Evf ~ ^ K^^^v) 



' 2xE^ ■ 2E^ ^-^^^ • 
We now rescale the Lagrange multiplier iVoid — NncwG{E^y / E'^ , the rescaled Hamiltonian constraint is, 

E[ = i?vac + 2G iJmatt (7) 

where 

iJvac = + =dH,{x)/dx, (8) 

„ _ x^ifT^xK^P^ 

^matt - 2(£;^)2 + 2{E'P)^ EV ' ^ ' 

This form of the Hamiltonian constraint allows an easy identification of the required boundary term if one assumes 
asymptotically flat conditions. The total Hamiltonian is given by, 

Ht= [ dxN{x){H,,,{x) + 2GH^,tt{x))+HB (10) 
Jo 

where N{x) is the rescaled lapse A^ncw and Hb is the boundary term at the asymptotic region x^ . Integrating by 
parts we get 

Ht = - dx^^^(^H4x) + 2G dyHmMy^j +Nix+)(-2GM + 2G dyH,,,^tt{y)\ + Hb 



r dx^^^ (n^ix) -2G [ dyH^^ttiy) + 2Gm] - 2GMt. 
Jo dx \ Jx / 



(11) 
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The boundary term Hb = —2GMt has been introduced in order to ensure that M is a constant and r the proper 
time in the asymptotic region. This is the standard boundary term in the spherically symmetric case. M is the space 

time mass while the Schwarzschild radius is given by Rs — 2G{M — dyH-caa.tt{v)))- In the case of a space time 
with a black hole the radial coordinate is given by i? = a; + Rs- M is a Dirac observable. In the case of weak fields 
therefore, so is the integral from to oo of iJmatt that we shall cal Hm- Even in presence of black holes Hm is an 
observable if the black hole is isolated. We will treat Hm as an energy in order to define the vacuum and the excited 
states of the theory in the case of interest in this paper, weak fields without the presence of black holes. 



III. QUANTIZATION OF THE MATTER FIELD ON A FIXED FLAT BACKGROUND 

Since we wish to understand in which way loop quantum gravity recovers results from ordinary quantum field theory 
in curved spacetime, we would like to outline some of those results for later comparison. If the space-time is flat it 
is convenient to fix the gauge K^p = to obtain explicitly the background metric in the usual spherical coordinates. 
In this case one solves i?vac — one gets that E"^ = x. Solving the evolution equation yields the Lagrange multiplier 
and one recovers the full flat space-time metric. The matter portion of the Hamiltonian constraint becomes, 

^matt-^ + • (l^j 

The evolution equation obtained from this Hamiltonian corresponds to spherical waves, 



2— =0. (13) 

X 



This can be solved by separation of variables. 



(C(tj)exp(-itjt)-f C(tj)exp(iwt))sin(cja;) 
0(x,t)=/ duj^ = '- , (14) 



which corresponds to spherical waves that are regular at the origin. From Hamilton's equation we can get an expression 
for P^, 

f°° {—iC{uj)ueyi\){—iut) + iC{u)ijjeyi]){iijot))xsin{uLix) 
P^(x,t) = / duj^ = . (15) 



From the standard commutation relations, [4>{x,t), P^{y,t)] = iS{x — y), one gets the [C{uj),C{ui')] = 5{uj — lu'). 
One can proceed to define a vacuum state |0) as the state that is annihilated by C. If one evaluates the expectation 
value of -ffmatt on the vacuum state one finds that it has an ultraviolet divergence. The usual resolution of this 
problem is to introduce a cutoff. It should be noted that when one treats this problem in loop quantum gravity this 
type of divergence does not appear because the well defined objects are holonomies associated to finite paths. In our 
treatment this aspect is lost since we have gauge fixed the radial variable which therefore becomes a c-number. As 
we usually proceed when we use the uniform discretization technique, we regularize the expression by placing it on a 
lattice. We will discuss later on the issue of taking the lattice spacing to zero. 

We will assume that the radial direction is bounded with a spatial extent L and consists of discrete points Xi 
separated by a coordinate distance e, and in particular we take Xi as e times an integer. We reinterpret the integrals 
as sums, Dirac deltas as Kronecker deltas, functional derivatives as partial derivatives, and partial derivatives in the 
radial directions as finite differences. Specifically Q 



/ 



dx ^ (16) 

X 

5{x-y) - ^ (17) 



5 1 d 



S<j){x) e 



(18) 



<^(x)' - (19) 
icof ^ i:.(2-2cos(.a;0) ^^^^ 
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If the spatial direction is discrete, the associated momentum space is bounded with extent 27r/e. To the first nontrivial 
order in in epsilon, all formulae involving momenta uj are unchanged except that momentum integrals are now sums 
over a momentum space of finite extent. 

The expectation value of i?matt can be computed replacing the quantum version of the expressions given above for 
t) and P^tix, t) in -ffmatt- Computing the expectation value on the vacuum state one is only left with contributions 

proportional to CC . On the lattice the result may be approximated in the limit of large L by the integral, 

(0|^„.au(.)|0) = r'^rfc. "'^'-^^"'^"^^"^"^"'^^"''^""^ - (21) 

Jo ZX TTLO 

The integral can be computed in closed form in terms of integral cosine functions. It is more useful to give an 
approximation for its value as an expansion in e, 

, - , , TT sin^(27ra;/e) ln(x/e) n^ , . 

m^.u{-m + + (22) 

The leading order in the energy density expansion is 7r/e^ which has the correct dimensions for an energy density in 
one spatial dimension, since we are only considering the radial mode of the scalar field. 

As in four dimensions, the energy of the vacuum gives rise to a cosmological constant if one allows the field to 
back-react on gravity. The nature of this constant is different, however in two dimensions |9(|- First of all, notice that 
if one had started from four dimensional gravity with a cosmological constant and imposed spherical symmetry, one 
can view the model as a 1 + 1 dimensional theory with a dilaton with a mass given by the four dimensional cosmological 
constant. That is, it does not produce a term that behaves like a cosmological constant in 1 + 1 dimensions. The 
vacuum energy, by contrast produces a constant term in the Hamiltonian constraint. Second, notice that even in 
vacuum ffvac already has a constant term in it. So the energy of the vacuum essentially operates as a rescaling 
of that constant term, which in turn can be absorbed by a rescaling of the radial coordinate. In four dimensions, 
if one chooses a Planck scale cutoff it implies that the radius of curvature of space-time becomes of the order of 
Planck length, which is clearly unphysical. In spherical symmetry the presence of the constant can be reabsorbed 
in a redefinition of the coordinates. This redefinition however, has consequences when one wishes to reinterpret the 
model as an approximation to a four dimensional space-time. The redefinition of the radial coordinate implies that 
the spheres do not have AttR^ area anymore. The four dimensional universe modeled contains a topological defect, a 
"global texture" [13] ■ Notice that this immediately precludes taking the lattice spacing to zero, since already when 
the lattice spacing is of the order of ^pianck one will have a solid deficit angle that exceeds 47r and does not allow to 
interpret the model as a four dimensional space-time. 

There are two avenues to handle the situation: either one rescales the radial variable and accepts that the model 
approximates four dimensional space-times with (large) topological defects, or one can modify the two dimensional 
model by adding a constant to the Hamiltonian constraint (a cosmological constant in 1 -I- 1 dimensional gravity). 
Such a model will not stem from a dimensional reduction of four dimensional gravity, but upon quantization will turn 
out to approximate four dimensional spherical gravity around a flat background without a topological defect. 

We will take the first point of view and write the Hamiltonian constraint as, H = i/vac + GTJmatt, where 

= (^~x{l~2A)~xKl + j^^ , (23) 
Mr_ xK^_ 

-Hmatt - ^^^-j2 + (^¥')2 ^ ^ ^ 

where A = ^Pvac and Pvac is the vacuum energy density. We choose h = c = 1 units. This rewriting of the constraint 
has the property that the expectation value of i/matt will be zero in the vacuum. 



IV. FULL QUANTIZATION OF THE MODEL 



We would like to write the master constraint based on the Hamiltonian constraint of the model we introduced in 
the last section. Although the discrete Hamiltonian constraint fails to close a first class algebra, we have showed 
in pj| that with the uniform discretization technique one can consistently treat the problem by minimizing the 
resulting master constraint. To write the master constraint at a quantum level we will polymerize the expression of 
the gravitational part of the constraint. We will not use a polymer representation in the scalar sector for simplicity 
and because we want to make contact with the usual treatments based on a Fock quantization. It is known that the 
Fock quantization for fields can be recovered from the polymer quantization [l^ . [l3| . 
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A. Variational technique to study the expectation value of the mcister constraint 



Here we will introduce a variational technique to minimize the master constraint. The technique is general, it 
is not restricted to the model we study in this paper. Wc start by considering a fiducial Hilbert space Tiaux in 
which the master constraint is a well defined self-adjoint operator. We will then use a variational technique to find 
approximations to the minimum value of the expectation value of the master constraint within this space. In many 
cases of interest, the minimum expectation value will not be zero, but will be small (the master constraint has units of 
action squared, so normally one would require it to be much smaller than , in order to have a good approximation 
of the physical space, in our units that translates into much smaller than one). As we will see in the examples, 
the resulting quantum theory will therefore not reproduce exactly the symmetries of the continuum theory but it 
will approximate them, even at the quantum level. We will see that if zero is in the spectrum of the operator the 
corresponding eigcnstates in many cases will be distributional with respect to the fiducial space we are considering. 

To implement the variational method, we consider trial states in Haux that are Gaussians centered around the 
classical solution of the model of interest in phase space. That means that as functions of Haux these will generically 
be Gaussians times phase factors such that the resulting state is centered around the classical solution in both 
configuration variables and momenta. The states are parameterized by the values of the standard deviations of the 
Gaussians in either configuration or momentum space. A caveat is that in gauge theories one may choose to work 
with a classical solution that is not in a completely determined gauge. Such a solution will be a trajectory in phase 
space. Such a trajectory will determine some of the canonical variables as functions of others, which will remain free. 
In that case one has to allow such variables to be free in the trial solution by considering Gaussians centered around 
a value that is a free parameter. If one chooses to work with a classical solution in a completely specified gauge 
one just considers Gaussians around the point in phase space represented by the classical solution of interest and 
extrcmizes the expectation value of the master constraint with respect to the standard deviations of the Gaussians. 
It can happen that the extremum occurs as a limit in the parameter space in which case the resulting state does not 
belong in Waux but in its dual (after a suitable rescaling, it becomes a distribution). 

Before attacking the problem of interest, it is useful to see the technique we just described in action in a couple 
of simple examples. The first example we choose is a system with two degrees of freedom qi,pi and q2,P2, and two 
constraints pi = and p2 = 0. The total Hamiltonian for the system is Ht = Nipi + N2P2 with A^i^2 Lagrange 
multipliers. The states annihilated by the constraints are trivial and given by the distribution 5{pi)5{p2). We fix a 
gauge qi — q2 = 0. Fixing the gauge is not needed in a simple model like this, but may be a necessity to simplify 
things in more complicated models. So we will choose a gauge fixing here to show that in the end the process loses 
all information about the gauge fixing and recovers the correct physical state. This requires fixing the Lagrange 
multipliers so there is only one (N) left and the total Hamiltonian becomes Ht = N{pi +P2)- The conjugate variable 
to the gauge fixing, pi — p2 is strongly zero. We start with a two parameter family of states in Waux choosing as 
configuration variables qi — 52 and pi + P2, 

V..,/, = ^^=exp (-^^i^) oxp (-^^i^) exp(i/3(pi +P2)) , (25) 

with P an arbitrary parameter associated with the fact that the variable qi + 32 is pure gauge. One could choose 
to work in a completely gauge fixed solution in which qi + q2 is zero, in that case there is no need to introduce the 
parameter f3. The choice of this family of states is based on the fact that they describe wave-packets centered around 
the classical solutions of the constraints, qi—q2 = 0, pi — P2 = and pi-\-p2 =0. We now define the master constraint 
M. = p\+p\ and act on this space of states. The expectation value is, 

(Va±,/3|H|V'.±,;3) = ^ + (26) 



where \/{cr±) are the standard deviations of the Gaussians, a± taken to be positive. One therefore sees that the 
expectation value cannot be zero for any finite value of the cr's. However, if one takes cr_ = and a+ = 2e^ then in 
the hmit e ^ 0, < H >= O(e^). The states \ip^) become, 

(91 - 92, Pi +P2|V'e) = ^ exp (gi - q2f e^) exp i - ^^^^J"^^ j exp {i(3 (pi +P2)) , (27) 



And their Fourier transform 



(pi -P2,Pi+P2\A) = ^:^^*^^P y ^^^^J^^ j y ^^^tj^^ j exp(i/?(pi +P2)) , (28) 
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These states are normalized in Tiaux but they vanish (in the sense of distributions) in the hmit e ^ 0. They need 
to be rescaled in order to end up with well defined distribution on some suitable subspace of Tiaux- 
So the physical states would be 

{Pi -P2,Pi +P2|V')ph = lini — ^(pi -P2,pi +P2|i/'e) = 2(5(pi +p2)5{pi -P2) = S{pi)6{p2) (29) 



Notice that the parameter (3 is free at the end of the process since it corresponds to the value of a variable that is 
pure gauge in this model. 

There is an additional element that the above example does not capture and we would like to discuss. When we 
apply this technique in situations of interest, we will be discretizing the theories we analyze. Usually, discretization 
turns first class constraints into second class ones. The uniform discretization procedure tells us that we do not need 
to concern ourselves with the second class nature of the constraints (for a discussion see [HI). We can still consider 
the master constraint and seek the minimization of its eigenvalues, but the presence of second class constraints in the 
discrete theory usually implies that the minimum eigenvalue of the master constraint will not be zero. The best one 
can hope for is that it will be small and the resulting quantum theory will approximate the symmetries of the theory 
one started with. This is a point of view that has been held as natural for some time in the context of quantum 
gravity, where one expects that some level of fundamental discreteness will emerge. We would like to illustrate this 
with a modification of the previous example. Instead of taking pi — and p2 — as the constraints we will take 
Pi + aq2 = and P2 = with a a small parameter (in realistic theories the small parameter is related to the lattice 
spacing in the discretization). We will still take the same set of tpc7±,i3 ^ts before, that is, for the trial solution we have 
chosen Gaussians centered around classical solutions of the gauge theory where the anomalous term vanishes. We do 
this because one usually knows solutions to the continuum theory one wishes to approximate (e.g. flat space or the 
Schwarzschild solution in the case of gravity) whereas the discrete theories have complicated solutions that usually 
cannot be treated in analytic form. The master constraint now becomes, 

M = pl+pl + 2apiq2 + a'^ql (30) 

and using the same ansatz (pS)) for the states one finds that 

We would like to identify a limit in the variables a± such that this quantity vanishes. As was to be expected, this is 
not possible. We can attempt to find values of the parameters a± and f3 that minimize this expression. The result is 
(3 = and (t+ = v^a and ct_ — which yields (■i/'min|IHI|V'min) = V^a. The state is, 

{puP2\i^n.n) = exp (.(eI^PlM] W ^. (32) 
\ a / y an 

It is interesting to compare this state and the corresponding expectation value of H obtained from our variational 
technique with the exact minimum of this model. A naive analysis would tell us that the minimum corresponds to 
an exact eigenstate with zero eigenvalue for H. However, that solution is not well behaved. It is known that one 
can find solutions of the master constraint that do not solve the constraints if one does not impose regularity in 
the solutions found The master constraint is an operator in the Hilbert space and one can analyze its spectral 
resolution. The spurious solutions do not belong in the spectral resolution of the master constraint. In this case 
one can solve exactly the eigenvalue problem M\^) = E\^p). The solutions with minimum eigenvalue are of the form 
<5(Pi)V'o(P2) where 4'o{P2) is the fundamental state of the Hamiltonian of a harmonic oscillator in the momentum 
representation. The minimum eigenvalue for such exact solution is a (compare with the variational one in which the 
eigenvalue was slightly higher \/2a). It is also interesting to note that if instead of choosing the gauge (71 — (72 = we 
had chosen qi = and proceeded with the variational technique, one obtains the exact state directly. This illustrates 
that the method approximates well the state of interest in situations where zero is not in the kernel of the master 
constraint. The solution that minimizes the master constraint admits a very simple interpretation that shows that the 
uniform discretization of the theory with the anomalous term a small but non- vanishing, approximately reproduces 
the invariances of the theory with first class constraints pi ~ P2 ~ 0. In fact qi and 52 are gauge variables and the 
physical space is independent of these variables. The physical state is constant in qi and 92- For a small but non 
vanishing alpha the physical states are independent of qi and weakly dependent on (72 . A final comment is that in this 
case the parameter /3, which was not determined in the case with first class constraints, gets determined here. That 
is, in the case where /? was associated with an exact gauge symmetry, the minimization of the master constraint was 
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insensitive to the value of (3. In the case where the constraints are second class and we do not get zero as minimum of 
the master constraint there is some dependence on /3, but it is weak, since the term in the master constraint is 0^0^ 
and a is small (in the quantum state one has approximately 5{p)exp{ipl3)). The theory where one does not exactly 
annihilate the master constraint only has approximate gauge symmetries and therefore has slightly "preferred" gauges 
from the point of view of minimizing the master constraint. 



B. The discrete master constraint 



Let us now consider the complete Hamiltonian constraint. We wish to discretize it and to polymerize the gravita- 
tional variables. The Hamiltonian gets rescaled in the discretization H{xi) H{i)/e,. We also rescale the expression 
multiplying the continuum Hamiltonian constraint times G. The resulting discrete expression is, 

Hi^) = -(1 - 2A). - + ,f^'ipKf + ^)) + + P^P^ " (33) 



G 



xpm? ^ x{^)-m + i)-m' ^ sin(px,(z)) _ 



We need to construct the master constraint. Since the Hamiltonian is a density of weight one, we define the master 
constraint associated with the Hamiltonian constraint in the full theory as, 

H = i 1 dx^^ip, (34) 

or, in terms of the variables of the model, up to a constant factor, 

H = 1 /dx J^^€p, (35) 



and in the discretized theory H'^ = with 



H(.)4^£i3— (36) 
2 y^E^E^ii) 



The constant Ip must be introduced so that H is dimensionless with h — c — 1, one could use vG instead of it. It 
is convenient to rescale the Hamiltonian constraint by E'f /{E'-^)' . This does not change the density weight. If one 
does not rescale things it turns out H is proportional to 1/E'^ . In the polymer representation this implies that the 
vacuum is the "zero loop" state, which is degenerate (it corresponds to zero volume space-times). To eliminate this 
unphysical possibility one exploits the fact that the Hamiltonian constraint is defined up to a factor given by a scalar 
function of the canonical variables without changing the first class nature of the classical constraint algebra. The 
rescaling factor in the discrete theory after the gauge fixing is E'p {i) / {2x{i)e) . So ([33]) has to be multiplied times 
that factor when constructing the master constraint ()36p . 

Let us focus on the matter portion of the Hamiltonian, we will write it as, 

rr l-\ ^^Ittji) , -^mitt(»)sin(P-^y(»)) „(3) ,.^ ,o-7^ 

HrnM^) = j^^^ + H^,,M. (37) 

The master constraint can be written as, 

2 



Cllit) (h^U^)) +C22(*)(i/,SttW) (38) 
~C,2m^Lm^2tt{i) + Ci3(*)i?i'itt(*)i?i'itt(0 + C23(*)i?i'itt(0i?mitt(0 + COo(*) 



where, 



^matt 



{€) = (e(F^(^))Vex(^)4(0(^ + l)-0(^))2)4 (39) 
H^2tS) = {~2x{z)m + l)-m)P^{i))ep (40) 
^?i'itt(0 = 2pvace4- (41) 
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To economize spaee, we will not give the classical expressions for the coefficients, since they can be readily obtained 
from the quantum expressions. 

In order to quantize the master constraint we need to choose a factor ordering. The expression of the master 
constraint is a sum of symmetric operators consisting of polynomials in i?"^ and sin(pK^), P'* and cf). We choose a 
factor ordering with the factors of E'^ are distributed symmetrically to the right and the left of the factors of sm{pK^p). 
For the factors f"^ and ^ we follow a similar strategy, putting the P*" symmetrically to the left and to the right of 
(f)'s. The coefficients in the above expression of the master constraint with this factor ordering are. 



Cll(«) 
Cl2(i) 
Cl3(i) 

C22(«) 
C33(i) 

ci(i) 



C2{l) 



4a;(i)2ei;v(i)4 
1 1 

2a;(i)2pe E^[if/2 
1 1 

1 / 1 



sm{pK^p{i)) 



iv(i)3/2' 



1 



8x(i)2p2e yEfiif Ef{i) 

sm{pK^{i)) 



- cos{2pK^{i)) 



E^ii) 



1 



4x(z)2e' 

x{i)e 



1 



1 cos{2pK^{i)) 1 



-2e(l-2A) 
1 



x{i + l) cos{2pK^{i)) x{i + 1) \ 1 



+ 



ex{i + If 



E'P{i) 



n2 



Ev{i) 4x(i)p'^e Ev{i) 4:x{i)p^eEv{i)2 2x{i) 
^ (1 - 2A) + .lit!]. {oosi2pK^{z + 1)) - 1) + 



(^E'P(i)Ef{i + l)^ 
3 exii + l)^ 



2px{i) 
1 



sm{pK^{i)) 



1 



1 



^x{i)e 



Ev{i) 



: s\n{ZpK^{i)) 



8p^xii)e 2pa;(i)2£;v(i + l)2 
1 



C3(i) 



coo(i) = 



E^{i) 
x{i)e 1 

2p E'^(i)3/2 

^ (1 - 2A) + (1 - cos{2pK^{i + 1))) 



sin(/9ii'<^(i)) 



Ev{ifn 



2x{i)-^ 
x{i)e 



Ax{i)'^ep'^ 
ex{i + If 



4a;(i)ep2 



(1 - cos{2pK^{i))) 



2E^{if 2x{ifEv{i + lY 

(3 - 4.cos{2pK^{i)) + cos(4pi4:^(i))) 



32ep' 



n2 ' jfe^2§(l-co«(2p^v'(^ + l))) 



4x(i)2 4x(«)2p' 
x{i + 1) 



8a;(i)ep4 



32f/94.-r(z)2 
Aa;(i + 1) 



(1 - cos(2pis:<^(i)) - cos(2pi4:^(i + 1)) + cos{2pK^{i)) cos{2pK^{i + 1))) 
(3 + cos{ipK^{i + 1)) - 4cos(2pii'y(i + 1))) 



2x(ifp^ (l-cos(2pi^^(i + l))) 

^ r(l-2A)(l-cos(2pX^(i))) 



4a;(i)p2 



eA 



(1-A) 



(42) 
(43) 

(44) 

(45) 
(46) 

(47) 
(48) 



(49) 



(50) 
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+ , ' cos{2pKJi))-. 

-44^ I ^(0 J—cos{2pK^{i))J—\ 

V \ \Ev{i)-i Ev{i) ^ ^^"Ev{i)J 

-x(i + 1) I J— cos(2pKJi + 1))^^^ I I 

^ ^ \Ev{i)^ Ev{i) ^ " E^ii) j j 

ex{i + l)'^ ( 1 1 /o /■ 1 \ a;(i)e2 

\ n ^ cofi{2pKJi + 1))^^ H .^-^^ (1 - 2A) 

4a:(z)2p2 \^£;v(i + i)2 + ' Ev'{i + l)) 2Ev{iy^ 

and it should be noted that the coefficients commute with H^^l^^, H^l^^ and -ff^^tt there are no ordering issues 
with them. 



C. Construction of the trial states 



Since we are interested in the vacuum solution, that classically corresponds to vanishing scalar fields, we will 
therefore ignore i?matt ((24|) and only consider the gravitational part ((23l) in order to construct the classical solution 
used to build the ansatz states for the variational technique. 



^^(l-2A)-:rif2 + __J . (52) 

As we discussed in subsection A we will choose a definite gauge to work in. Our choice is A'<^ = 0, and this implies 
E'^ = xl\j\ — 2A. As we claimed before, the presence of the cosmological constant rescales the radial variable (recall 
that without the constant the solution was E'-^ = a;). The resulting four dimensional space-time will be locally flat 
with a solid deficit angle and described in spherical coordinates. 

We construct a polymer representation. As we did in previous papers 5] one sets up a lattice of points j — . . . N 
in the radial direction and writes a "point holonomy" for the variable at each lattice site, 

Tp = exp j = {K^lfl). (53) 

In this expression the quantities pi play the role of the "loop" in this one dimensional context. They also are 
proportional to the eigenvalues of the triad operator E''^{i). The quantum state we will choose for the variational 
method will be centered around the classical solution and therefore we will choose to have the variable pi centered at 
the classical value of £"^(i) = exi(i) = ex{i)/\/l — 2A, 



on this state (£"^(i)) = tx\{i) and (K^p(i)) = 0. Notice that this type of ansatz in general will be too restrictive: 
we have ignored possible correlations among neighboring points by assuming a Gaussian at each point. This could 
potentially be problematic when studying excited states and computing propagators. We will not attack those 
problems in this paper so we will continue with the restrictive ansatz for the moment being. 

We will now compute the expectation value of the matter portion of the Hamiltonian constraint on the above 
state. The result will be an operator acting on the matter fields. We will then construct the vacuum for the resulting 
operator. What we are doing is to construct a quantum field theory living on the geometry given by the expectation 
values of the triad and extrinsic curvature on the above state. We proceed in this way for expediency since this is our 
first approach to the problem. In the future we plan to revisit the problem treating all the variables in a polymerized 
representation, both gravitational and material ones, with the variational technique. Preliminary results indicate that 
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such an approach is viable. For the matter field one would start by considering a coherent state centered at zero 
values for the field and then will obtain the vacuum as a limit. This would yield valuable insights into the relation of 
the usual Fock quantization with the loop quantum gravity techniques, especially when one gets to discuss physical 
elements like the propagators of fields. 

In order to take the expectation value of the matter portion of the Hamiltonian constraint, ([37)1 on the state ([54]) 
we need to realize two quantum operators. The first one is. 



where we have considered the action on one of the factors of ([54]) . To derive this expression we consider 

/ ,N -3/2 . . .X -3/2 . ,N -3/2 

lEj E iEj and use the realization of iEj that was discussed in the context of loop quantum cos- 
mology in [i5[. The reason we can use the loop quantum cosmology results is that our Hilbert space is a direct 
product of loop quantum cosmology Hilbert spaces each at one of the lattice sites in the radial direction. With the 
above result one can compute the expectation value. 



1 ^ 1-2A , 5 ep^{l-2Afp^ ^ 3 aep^{l-2Af 



The calculation is done by integrating in fl and the result is lengthy, here we just show it in the approximation 
e > i-p. The first term is the classical value, the others are quantum corrections, the first one comes from the 
polymerization, the second from fluctuations in jl. The second operator we need is the one arising in the second term 
of the Hamiltonian, 

/ / I 1 sm{pk^{i)) 1 

(V'ffI , = I yg) = 0- (57) 



To quickly see why this is zero keep in mind that the state is a Gaussian centered at ~ and the sine is an odd 
function. With these results the expectation value of the Hamiltonian (the "effective Hamiltonian" ) is, 

/ff,,, = (V^.H^matt(2:,OI V-.^) - + '- ~ Pvac (58) 

In this equation we have pursued the unusual approach of taking the continuum limit in the terms that involve 
derivatives and the terms that involve the momenta of the scalar field. This simplifies calculations since we will be 
dealing with differential equations rather than difference equations. The idea is that the solutions to the differential 
equations, suitably discretized, will be a good approximation (at least to 0{e) corrections) to the solutions of the 
difference equations. In the above expression the quantity g{x) is given by, 



5 £pV(1-2A) 3g(x)4(l-2A) 
16 a;2e2 ~ 8 

From the effective Hamiltonian we get the "wave equation" for the fields living on the curved semiclassical back- 
ground, 

2 d<P{x,t) 2 d<j>{x,t)dg{x) d'cj,{x,t) 1 gjxf d^cj,{x,t) _ 

X dx g{x) dx dx dx^ 4 (1 - 2A)2 dt^ ' ^ ' 

Since the background is time-independent, positive and negative frequency modes can be introduced by going to 
Fourier space in t. The resulting equation can be cast in Sturm-Liouville form as, 

^2 

{2B{x)(t)'{x,u;))' + = (61) 

where 

Af.) - ^^-^^-^^ (62) 
^^"^^ " 1-2A 8 62 4 e2 ' (62) 

Bix)=xHl-2A) + l'4^ + ',^-^ (63) 
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The solution to this Sturm-LiouviUe problem is 

<l>{x,w) = -sin ( — - ) (64) 



X V2(1-2A) 



1 

'3^ 



X uo COS ^1 (^2;) — —u cos (^~^ J ^ a; w sm (-^H^^;) + sm y~j 



4e2 



— + 3. 



and this solution neglects terms with higher powers than £p/{€x)^. Where Si(a::) = f^dtsin{t)/t, and Ci(a;) = 
7 + ln(x) + dt{cos{t) — l)/t are the sine integral and cosine integral functions respectively and Euler's Gamma 
is given by 7 = 0.5772156649. The first term in the bracket in ([64|) corresponds to the standard spherical mode 
decomposition in (locally) flat space-time. The next parenthesis includes two terms that are corrections, the first due 
to polymerization and the next, involving cr is a quantum correction. These terms would not be present in a treatment 
of quantum field theory on a classical space-time. Using the Hamilton equations we can compute , 

^ 2V^(1-2A) dt ^^^^ 

and use it to compute the effective Hamiltonian ((SS]) . 



Hf,,, = (1 - 2A) ckuLodiu;)CiLo). (66) 

Jo 

To obtain this expression we note that the solution (j64p can be written as 0(a;, t) = duju{x, uj)h{uj, t) where hiuj, t) 
is the last parenthesis in ([M]) . Notice that we have introduced a lattice cutoff for the frequency 27r/e. Then one 
uses the lattice version of the closure relation dLJu{x, uj)u{x' , ui) — 2S{x — x')/A{x) and the orthogonality relation 
dxA(x)u{x,uj)u{x,uj')/2 = S{uj — lo'). 
We have therefore concluded the computation of the state that we will use as a trial in the variational method. It 
will be given by a direct product of the vacuum of the matter part of the Hamiltonian (|66p and the Gaussian (|54p on 
the gravitational variables. 

\^f-^) = \rp,) ^ |0) (67) 

The parameters a will be varied to minimize the master constraint. Notice that the state is a direct product because 
we are considering the vacuum. If we were to consider excitations then there might be entanglement between the 
matter and gravitational variables [T6| . 

D. Minimizing the master constraint 



The realization of the master constraint psp as a quantum operator depends on the realization of six key operators. 
We proceed to present their expectation values here. We start by the operators involving the cosine of K^p, 

COS (2pk^{i)) l^r^i) = exp (-^) , (68) 

(V^r'^'l COS [^pKii)) W'^') = exp (y-^^ ■ (69) 

We then consider the powers of the inverse of E''^ . We already computed the expectation value of the square in 
Here we list the other needed powers, 



1 t^i^i (1-2A)2 5^p^(l-2A)>2 5(j£p4(i_2A) 



(71) 



(^^"''It ^ [pKki)) ^1^^""') = 0, (72) 
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(V4""1^^ sin (pkA^)) ^^|V>r^') = 0, (73) 



(^r^^l^^ sin (ipk^ii)) — ^IV^r^^i) = 0. (74) 



With these results we can proceed to compute the expectation value of the master constraint on the gravitational 
state. The result will be an operator acting on the matter part. The calculation of the expectation values of the 
coefficients Ci and (|42 |) - (f5T|l is straightforward, but lengthy. We will not list the results here. What is more 
challenging is the computation of the expectation value of the matter part of the expansion of ([551) . It helps that 
some of the coefficients vanish. The non-vanishing contributions are, 



2 



(V5|H(i)|^5) = 



(75) 



We now need to compute the expectation value of this operator on the matter vacuum. To do this we again 
use the procedure of going to the continuum limit in the matter terms involving derivatives and momenta and 
integrating in the frequencies with an ultraviolet cutoff. Let us start with H^l^^{i). The continuum limit expression 
is H^l^^{x, t) = £p ^(P'^(x, t))^ + X** (0'(x, <))^^ . We now substitute P'^ and 4> by their mode decomposition and one 

gets a quadratic expression in the C"s and u's. The expectation value only gets contributions from the CC terms. 
The result is. 



and substituting u{lu,x) and A{x) we obtain, 

. > X . ^9 / 1 cos^(^) 7rsin(2^)cos(^)\ 

m^Uxm = /^(i-2A)A(..)^(^ + ^-^ 4^7^) 

, I'p /tt^x^ ln(2) a;^cos(^)sin(^) 5 . ^.ttx^ ^ 1 . .ttxA 

where Cin(x) = 7 + Incc — Ci(a;). One can get a more manageable expression, which we will use in the rest of the 
paper by ignoring corrections of ip and neglecting the highly oscillating terms that involve sin(7ra;/e) or cosines and 
the integral cosines. The result is, 

m^^Uxm = ^(Y^ (-2 + ^ + ln(2) + 7 + ln(^)) , (78) 
and the dominant term is Tr^z^/e^. Reverting to the discrete theory, it reads, 

r(l) f,Mn\ _ '''k^ ( o , T^^^iif . . .^2: 



The procedure to compute the expectation value of the other terms in ((75|) is exactly the same, but the size of the 
expressions involved is quite large. We will not display them here for reasons of space. 
The result for the expectation value of the integrand of the master constraint is. 



ex^ 



ej ex\l- 2A) 96 e^xVo (1 - 2A)' 



1 Att^ 43 Att 



e 



SxVo (1 - 2A)' 128 e^x^ao (1 - 2A)' 



e(7-2 + ln(2^)). , ^ 2vr ,^fL,,,_ 



c4(l-2A)L2 a;2^4 ea;4(i_2A) \ej ex^L^ i2ex^{l~2A) 
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FIG. 1: The expectation value of the master constraint as a function of the lattice spacing. We see that the value of the 
master constraint is small unless one chooses lattice separations of order Planck length. The figure does not show it, but for 
separations of the order of 10~^^cm the master constraint is very small, of the order of 10^^" (we are using units in which h is 
one and therefore the master constraint is dimensionless) . 




1 TT^ 32e / /L 

+ T7: —r. + In 



x^L^ V e / 2:3^2 256 £33.4^^ (1 _ 2A)" cj^tx^iT 



1 (4x^6^0 + 4 xaoV + 4ao6^x-H7aoV) /I i f 2 + \J'^^''W\ 

4(1^ {x + e)alx^e^ ^ 1,^ ^ V ' ^ \~ ) ) ) 



O-OTT , ^ CToTT 6 (To 



{x + e)x^e^ {x + e)x'^L^ {x + e)x^n \e J x^eaoL^ ex^{l~2A)L^ 



-2 ^TT-^ITT^ ( 7 - 2 + In I — ) ) In I - ) + 3 ::f. + 6 In I - ) - 1/16 ^ 

We have assumed u — aoe^/£p with cto of order unity and we have neglected terms 0{£p). We have assumed a to be 
independent of x in order to simphfy the above expression, which otherwise becomes too large. Experiments we have 
carried out suggest that allowing variations in x leads to the same minimum value of a approximately independent 
of X. 

We would like to study the minimum of the master constraint as a function of (Jq for different choices of e/ip. Notice 
that we have assumed ctq to be of order one. One can change that by varying the ansatz for a including other powers 
e/£p different than 2. We have carried out such experiments. The results can be summarized as follows. In figure 
[T]we show the value of the master constraint as a function of e (in centimeters) and for cto — 10 and cr = <Joe^/£p. 
Varying ctq while keeping it of order one changes little the shape of the curve. We see that in the approximation 
studied the theory does not appear to have a continuum limit, but we see that the master constraint quickly drops 
to zero for lattice spacings larger than Planck scale. Although the figure suggests that the master constraint drops 
even further for larger lattice spacings, the approximation in which we have handled expressions (in which we have 
neglected higher powers of e/(p) is inadequate for large values of e and the master constraint very likely will increase 
its value for large values of e. So there exists a genuine preferred value of e that minimizes the master constraint. 
Even so, the approximation should be reliable up to values of e ~ 10~23cm and for such values the master constraint 
is of the order of 10"^'', so one sees that this is a regime where one approximates the continuum theory very well. 

We have explored other ranges of cr's (with different powers of e/^p). The observation is the following. For lower 
powers than three we get a curve that looks similar to the one shown in the figure, but that grows faster as one 
approaches smaller lattice spacings and therefore the minimum occurs farther away from the Planck scale. For powers 
higher than 10/3 one violates the approximation that £p /e is small and the expressions we derived are not valid. From 
these considerations and an analysis of the powers involved, we conclude that the minimum for the master constraint 
is achieved for a power of e/£p in cr close to two and e ~ 10^'^^p. 

An interesting speculation is that if the minimum of the master constraint happens in the range mentioned, the 
cosmological constant, which goes as A ~ ^p/e^ would not be of Planck scale but several orders of magnitude smaller. 

Another observation of interest is to note what would have happened if instead of choosing the state peaked around 
the fiat metric (with a topological defect) one would have chosen the "loop quantum gravity vacuum" , i.e. a state 
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with zero loops which corresponds to a degenerate metric \fJ.(i) = 0). Such a state annihilates the matter Hamiltonian 
in the loop representation and has zero volume. It would be disturbing if this state yielded a lower value for the 
master constraint than the state we constructed, since it would imply that degenerate geometries dominate. This is 
not the case, as can be easily seen. For such a state all expectation values ([55|) - ([7^ vanish. One can check that the 
expectation value of the master constraint is, 

m^l^^ (SI) 

That is, the result is very large. For e ~ £p it goes as L/£p, the size of the universe in Planck lengths. Therefore 
these degenerate states are heavily suppressed. 



V. DISCUSSION 



We have studied spherically symmetric gravity coupled to a spherically symmetric scalar field using loop quantum 
gravity techniques. The problem has a non-Lie algebra of constraints and we used the "uniform discretization" 
technique to treat the dynamics. We used a variational technique to minimize the discrete master constraint. With 
the trial states proposed, we were not able to reach a zero eigenvalue for the master constraint, that is, the theory 
does not seem to have a quantum continuum limit. The lowest eigenstate of the master constraint has the form of 
a direct product of a Fock vacuum for the scalar field and Gaussian states centered around flat space-time for the 
gravitational variables. Although the theory does not have a continuum limit, it approximates general relativity well 
for small values of the lattice separation, which in turn regularizes the cosmological constant. The lattice treatment 
we have performed diverges when one takes the continuum limit. The reader may wonder why loop quantum gravity 
has failed to act as the "natural regulator of matter quantum field theories" as claimed, for instance in [17]. The 
problem arises with the gauge fixing of the diffeomorphism constraint that we performed at the classical level. This 
leads us to variables that have the structure of a Bohr compactification in the "transverse" ip direction, but the 
variable in the radial direction is a c-number and therefore is not dynamical and has continuous character. There is 
no chance therefore that loop quantum gravity based on this gauge fixing could regulate the short distance behavior, 
which is responsible for the emergence of the cosmological constant. To tackle this issue one would have to allow 
both the diffeomorphism and Hamiltonian constraint to remain in the theory. The calculational complexity would 
increase importantly, since one will have to regulate the master constraint in such a way that the resulting states have 
remnants of diffeomorphism invariance in the discrete theory. This has been successfully accomplished with uniform 
discretizations in the Husain-Kuchaf model [11], but the complexity there was considerably reduced by the lack of a 
Hamiltonian constraint. It is worthwhile noticing that even if one allowed loop quantum gravity to regulate matter 
in the proposed way, the resulting cosmological constant is likely to be finite but still very large with respect to the 
current observed value. 

The present paper is a first exploration of a difficult problem, carried out with several assumptions and limitations 
that we have outlined in the text. Future work will include relaxing the assumption that one has a Fock vacuum for the 
scalar field and treating both the gravitational and scalar variables on the same footing with the variational technique 
for the master constraint. In this context it will be interesting to study the excited states of matter and study the 
modifications in dispersion relations for the matter fields due to the quantum geometry. This will definitely require 
considering trial states with correlations in the variational method, something we have not done here. One should also 
relax the gauge fixing of diffeomorphisms to see if the cosmological constant problem becomes better under control. 
Other future directions would be to consider solutions centered around non-flat geometries, for instance, including a 
black hole with the aim of studying if the scalar field states involve Hawking radiation. 
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